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1. Introduction 

Deep inelastic scattering (DIS) of electrons and positrons on protons at HERA has been cen¬ 
tral to the exploration of proton structure and quark-gluon interaction dynamics as described by 
perturbative Quantum Chromo Dynamics (pQCD). HERA was operated at a centre-of-mass en¬ 
ergy of up to ~ 320GeV. This enabled the two collaborations, HI and ZEUS, to explore a 
large phase space in Bjorken x, xbj, and negative four-momentum-transfer squared, Q^. Cross 
sections for neutral current (NC) interactions were published for 0.045 < Q^< 50000 GeV2 and 
6 • 10-^ < XBj < 0.65. 

HERA was operated in two phases: HERA I, from 1992 to 2000, and HERA II, from 2002 
to 2007. It was operated with an electron beam energy of Ee ~ 27.5 GeV. Eor most of HERA I 
and II, the proton beam energy was Ep = 920 GeV, resulting in the highest centre-of-mass energy 
of ^/s ~ 320GeV. During all of HERA running, the HI and ZEUS collaborations collected total 
integrated luminosities of approximately 500 pb^^ each, divided about equally between e+p and 
ep scattering. The data presented here is the final combination of HERA inclusive data based on 
all published HI and ZEUS measurements corrected to zero beam polarisation. This includes data 
taken with proton beam energies of Ep = 920, 820, 575 and 460 GeV corresponding to at ^/s ~ 320, 
300, 251 and 225 GeV. 

The combination was performed using the package HERAverager [|]] and the pQCD analysis 
using HERAEitter |^. The correlated systematic uncertainties and global normalisations were 
treated such that one coherent data set was obtained. The combination leads to a significantly 
reduced uncertainty compared to the orginal inputs and compared to the previous combination 
of HERA-I data, particularly in the electron sector, see Eig. ||. The combined data demonstrate 
electroweak unification beautifully and allow an extraction of xE^, see Eig. ^ 

Within the framework of pQCD, the proton is described by parton distribution functions 
(PDEs) which provide probabilities for a particle to scatter off partons, gluons or quarks, carrying 
the fraction x of the proton momentum. Perturbative QCD determines the evolution of the PDEs 
to any scale once they are provided at a starting scale. The name HERAPDE stands for a pQCD 
analysis, within the DGEAP formalism, to determine the PDEs at the starting scale by fitting the 
xb] and dependences of the combined HERA NC and CC DIS cross sections. The name HER- 
APDE2.0 refers to this analysis based on the newly combined inclusive DIS cross sections from 
all of HERA I and HERA II. The strength of the HERAPDE approach is that one coherent high- 
precision data set containing NC and CC cross sections is used as input. The newly combined data 
entering the HERAPDE2.0 analysis span four orders of magnitude in and XBj- The availability 
of precision NC and CC cross sections over such a large large phase space allows HERAPDE to 
be based on ep scattering data only and makes HERAPDE independent of any nuclear corrections. 
The difference between the NC e^p and e^p cross sections at high Q^, together with the high-Q^ 
CC data, constrain the valence quark distributions. The CC data also constrain the down sea-quark 
distribution in the proton without assuming isospin symmetry. The lower-Q^ NC data constrain the 
low-x sea-quark distributions. The precisely measured variations of the DIS cross sections in 
different bins of xbj constrain the gluon distribution. Measurement of cross sections at different 
beam energies constrains the longitudinal structure function, El, and thus provides independent 
information on the gluon distribution. 
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Figure 1: HERA combined NC e^p reduced cross sections as a function of for selected xej-bins com¬ 
pared to the individual HI and ZEUS data(top); and HERA combined NC e^p reduced cross sections com¬ 
pared to the to the HERA-I combination (bottom). 
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Figure 2: NC and CC e^p and e^p cross sections (top) The structure function at = lOOOGeV^ 
(bottom). The data are compared with the prediction from HERAPDF2.0 NLO. 
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The consistency of the input data allowed the determination of the experimental uncertain¬ 
ties on the HERAPDF2.0 parton distributions using rigorous statistical methods. Uncertainties 
resulting from model assumptions and from the choice of the parameterisation of the PDFs are 
considered separately. 

Both HI and ZEUS also published charm-production cross sections, which were combined and 
analysed previously, as well as jet-production cross sections. These data were included to obtain a 
variant HERAPDF2.0Jets. The inclusion of jet cross-sections made it possible to simultaneously 
determine the PDFs and the strong coupling constant o:, 5 (M|). Full details of the analysis are given 
in ref. [||]. 

2. HERAPDF2.0 and its variations 

Fig. H shows summary plots at /ij = lOGeV^ of the valence, total Sea and gluon PDFs for 
HERAPDF2.0 analysed at NFO and at NNFO, for the standard cut > 3.5GeV^. The experi¬ 
mental uncertainties are shown in red. Model uncertainties are shown in yellow. These are due 
to variation of the central choices for: the cut; the values of the pole-masses of the charm and 
beauty quarks; the fractional contribution and the shape of the strange-PDF. HERA data on charm 
and beauty production determine the central choices of the heavy quark masses and their model 
variations. Parametrization uncertainties are shown in green. These are due to: variation of the 
starting scale; addition of extra parameters. The central choice of parametrization is determined as 
usual for HERAPDF analyses by saturation of the but the addition of extra parameters some¬ 
times results in close-by but distinct minima. Additionally seen on these figures is the result of an 
alternative gluon parametrisation HERAPDF2.0AG, for which the gluon must be positive definite 
for all above the starting scale. These PDFs are similar in to the standard ones at NFO but 
are disfavoured at NNFO. An FO set is also available using the alternative gluon parametrisation. 
It is shown compared to the NFO set in Fig. The standard fits use a value of the strong coupling 
constant, as{M^) = 0.118, at NFO and NNFO, and, = 0.130, at FO, but sets using a range 

of values from 0.110 to 0.130 are also available. 

A more extreme variation of the cut, > lOGeV^, has also been considered, resulting 
in the HERAPDF2.0HiQ2 PDFs, since it was observed that the x^ per degree of freedom of the 
fit decreases steadily until > lOGeV^. This is true for both NFO and NNFO fits and it is true 
independent of the heavy quark scheme used to analyse the data, see Fig. ||. In fact it depends 
mostly on the order to which Fi is evaluated. The fits do not favour the evaluation of Fi to 0{aj). 
It is also somewhat counter-intuitive that the x^ is not improved when going from NFO to NNFO 
within the same scheme, the fits do not favour the faster NNFO evolution. 

HERA kinematics are such that low is also low x. Thus HERAPDF2.0HiQ2 PDFs are used 
to assess any bias resulting from the inclusion of low-Q^, low-x data which might require analy¬ 
sis beyond the DGFAP formalism, such as: resummation of ln{l/x) terms, non-linear evolution 
equations and non-perturbative effects. Figs. show that there is no bias at high scale due to the 
inclusion of the lower data. 

Figs ^ compare the HERAPDF2.0NFO fit to the HERAPDFl.ONFO fit. One can see the 
reduction in the high-x uncertainties and the fact that the high-x Sea is now much less hard. Figs. || 
make the same comparison for the HERAPDF2.0NNFO and the HERAPDF1.5NNFO fit. Again 
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Figure 3: The parton distribution functions of HERAPDF2.0 NLO(top) and NNLO(bottom), xuy, xS = 
2x{U + D), xg, at fif = lOGeV^. The gluon and sea distributions are scaled down by a factor of 20. The 
experimental, model and parameterisation uncertainties are shown separately. The dotted lines represent 
HERAPDF2.0AG NLO with an alternative gluon parameterisation. 
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Figure 4: HERAPDF2.0 PDFs at LO and NFO are compared using experimental uncertaities only. The 
alternative gluon form of the parametrisation is used for both. 


the high-x uncertainty is reduced and in particular, the high-x gluon has a much reduced uncertainty 
band and its central value moves towards the lower end of the HERAPDF1.5 uncertainty band. 

Two sets of PDFs using a fixed flavour number scheme have been extracfed, as shown in Fig 
These differ from each ofher in fhree respecfs: fhe order af which Fl is evaluated 0{aj) (FF3A), 
0{as) (FF3B); whefher or nof runs wifh 3-llavours (FF3A) or wifh variable flavour (FF3B); and 
fhe use of pole masses (FF3A) or currenf masses (FF3B). 


Heavy flavour dafa from fhe charm combination has also been added fo fhe fif, buf if does nof 
make much difference once ifs consfraining effecf on fhe charm mass has been faken info accounf. 
Adding dafa on jef producfuion also does nof make much difference IF fhe value of a^Ml) is kepf 
fixed. However if as{M^) is free fhen jef dafa have a dramafic effecf in consfraining ifs value, 
see Fig. IF, where fhe profiles vs as{M^) are shown for fhe NFO and NNFO fifs fo inclusive 
dafa alone and fhe same profile is shown for fhe NFO fif including Jefs. (Note fhaf we cannof 
include jefs in an NNFO fif since jef production cross sections in DIS have nof been calculated fo 
NNFO). A simulfaneous til of fhe PDFs and fhe value of as{M^) can be made once fhe jef dafa 
are included resulting in fhe value, as{M^) = 0.1183 ±0.009(cx/7) ± 0.0005{model/param) ± 
0 . 0012 (/iad)j(QQQ 3 Q(scaZe), where ’had’ indicates exfra uncertainties due fo fhe hadronisafion of 
fhe jefs. The gluon PDF is sfrongly correlated fo fhe value of as{M^) and fhus, in a fif where 
as{M^) is free, fhe gluon uncertainly increases. However, provided lhal jef dafa are included in fhe 
fif fhis increase is nof dramatic, see Figs |^. 
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Figure 5: The dependence of /d-o.f on for HERAPDF2.0 fits: (top) using the RTOPT and FONLL 
schmes at NFO and NNFO; (bottom) using RTOPT, ACOT and FONFF-B schemes and fixed flavour number 
schemes at NFO. 


3. Conclusions 

The HI and ZEUS data on inclusive neutral and charged current cross sections have 
been combined into a data set with a total integrated luminosity of ~ \fb^^. This data set spans 
six orders of magnitude in both r and Q^. The combined cross sections were used as input to a 
pQCD analysis to extract the parton distribution functions HERAPDF2.0 at EO, NEO and NNEO. 
The effect of using various different heavy flavour schemes and different cuts on the data was 
investigated. All heavy flavour schemes show some sensitivity to the minimim cut, however 
the choice of this cut does not bias data at high scale significantly. Eor the standard fits the value 
of as{M^) is fixed, but a measurement of can be made if jet data are included in the fit, 

resulting in the value as{M^) = 0.1184 ±0.0016 at NEO, excluding scale uncertainties. The data 
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Figure 6: The parton distribution functions of HERAPDF2.0, xm,,, xdv, xS = 2x{U 
lOOOOGeV^ with = 3.5 GeV^ compared to the PDFs of HFRAPDFHiQ2 with Q 
for the NLO fits; bottom for the NNLO fits 
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+ D), xg, at jif = 
= lOGeV^. Top 
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Figure 7: The parton distribution functions of HERAPDF2.0 NLO, xm,,, xdy, xS = 2x(U + D), xg, at fif 
lOGeV^ compared to HERAPDFl.ONLO on log (top) and linear (bottom) scales. 
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Figure 8: The parton distribution functions of HERAPDF2.0 NNLO, xuy, xdv, xS = 2x{U + D), xg, at 
= lOGeV^ compared to HERAPDF1.5 NNLO on log (top) and linear (bottom) scales. 
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Figure 9: The parton distribution functions of HERAPDF2.0FF3A NFO and HFRAPDF2.0FF3B NFO, 
xuv,xdv,xS=2x(U +D),xg, at jif = lOGeV^. The experimental, model and parameterisation uncertainties 
are shown separately. 
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Figure 10: vs. for pQCD fits with different using data on (a) inclusive, charm and jet 

production at NLO, (b) inclusive ep scattering only at NLO, and (c) inclusive ep scattering only at NNLO. 


and the PDFs are available on https://www.desy.de/hlzeus/herapdf20. 
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Figure 11: The parton distribution functions of HERAPDF2.0Jets NLO, xuy, xdv, xS = 2x{U + D), xg, at 
Ilf = lOGeV^ with fixed as{M 2 ) = 0.118 (top) and free (bottom). The experimental, model and 

parameterisation uncertainties are shown separately. The hadronisation uncertainty is also included, but it is 
only visible for the fit with free as{M^). 
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